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ViO ' Abstract. We discuss the correlation between the evolution of the TeV emission and X-ray radiation observed in high-energy 

peaked BL Lac objects. We describe such a correlation by a simple power law FtcvCO f In the first part of this work 

we present correlations obtained for the activity of Mrk 501 observed in 1997 April and for the activity of Mrk 421 observed in 
2000 February. Our results obtained for Mrk 501 show that the index of the correlation (x) may strongly depend on the width 
and position of the spectral bands used for the comparison. The result of the correlation which we have obtained for Mrk 421 
is not informative. However, we discuss results of similar correlation obtained for this source by other authors. They report an 
] almost quadratic (x ~ 2) correlations observed between the evolution of the TeV and X-ray emission. In the second part of this 

■ paper we present a phenomenological model which describes the evolution of the synchrotron and inverse Compton emission 

of a simple spherical homogeneous source. Neglecting the radiative cooling of the particles we derive analytical expressions 
, that describe the evolution. Then we use a numerical code to investigate the impact of radiative cooling on the evolution. We 

' show that different forms of correlations can be obtained depending on the assumed evolution scenario and the spectral bands 

used for the calculation. However, the quadratic correlation observed during the decay phase of the flare observed in Mrk 421 
^ ^, on 2001 March 19 appears problematic for this basic modeling. The quadratic correlation can be explained only for specific 

I ' choices of the spectral bands used for the calculation. Therefore, looking for more robust solutions, we investigate the evolution 

O I of the emission generated by a cylindrical source. However this model does not provide robust solutions for the problem of a 

quadratic correlation. In principle the problem could be solved by the TeV emission generated by the self Compton scattering 
^ ' in the Thomson limit. However, we show that such a process requires unacceptably large values of the Doppler factor. Finally 

we briefly discuss the possible influence of the light travel time effect on our results. 
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1 . Introduction (hereafter IC) from the same electron population (although dif- 
ferent scenarios have been proposed, see e.g. Mannheim 1 9931 

The emission from BL Lacertae objects is dominated by the AharonianlOOO, Pohl & SchlickeiserlMI Mucke et aU2(m. 
intense relativistic boosted non-thermal continuum produced 
within a relativistic jet closely aligned with the line of sight, 

making these objects (together with Flat Spectrum Radio BL Lacs are further subdivided in Low-energy Peak BL 

Quasars, the other subgroup of the blazars family) the best lab- Lacs (LBL), exhibiting the synchrotron peak in the IR-optical 

oratories to study the physics of relativistic jets. region of the spectrum, and High-energy Peak BL Lacs (HBL), 

The overall emission, from radio to y-rays (extended in in which the synchrotron peak can he in the UV-X-ray band, 

some cases to the multi-TeV band), shows the presence of two The almost featureless optical spectra of HBL clearly indicate 

well-defined broad components, the first one peaking in the that the environment external to the jet is quite poor in soft 

optical-soft-X-ray bands, the second one in the GeV-TeV re- photons, suggesting that the high-energy emission is mainly 

gion. The low energy peak is attributed to synchrotron emis- produced through the Comptonization of the synchrotron ra- 

sion by relativistic electrons in the jet, while the second com- diation (Synchrotron Self-Compton emission, SSC). The few 

ponent is commonly believed to be Inverse Compton emission extragalactic TeV sources firmly detected so far (Mrk 421, 

Mrk 501, PKS 2344H-514, P KS 215 5-304, 1ES1959-H65, IBS 

Send offprint requests to: Krzysztof Katarzyhski 1426-1-428, e.g. Krawczynski|2004|for a recent review) belong 

kat@astro.uni.torun.pl to this class. 
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Since the discovery of the first BL Lac object emitting TeV 
radiation, Mrk 421, (Punch et al. 1992 ) TeV blazars have been 
the target of very intense observational and theoretical investi- 
gations. Indeed the possibility of observing the emission pro- 
duced by very high energy electrons (up to Lorentz factors of 
the order of 10^) coupled with observations in the X-ray band, 
where the synchrotron peak of these sources is usually located, 
offers a unique tool to probe the processes responsible for the 
cooling and the acceleration of relativistic particles. 

Studies conducted simultaneously in the X-ray and in the 
TeV bands are of particular importance, since in the simple 
SSC framework one expects that variations in X-rays and TeV 
should be closely correlated, being produced by electrons with 
similar energies. Assuming a typical magnetic field intensity 
B ~ 0.1 G and Doppler factor 6 ~ 10, photons with en- 
ergy E ~ 1 keV are emitted by electrons with Lorentz factor 
y ~ 10^. The same electrons will upscatter photons at energy 
E ~ ymc^ ~ 1 TeV (since the Lorentz factor is extremely large 
the scattering will occur in the Klein-Nishina regime even with 
optical target photons). In fact observations at X-ray and TeV 
energies (Buckley et al. 119961 Catanese et al. I1997> yielded 
significant evidence of correlated and simultaneous variability 
of the TeV and X-ray fluxes. During the X/TeV 1998 cam- 
paign on Mrk 421 a flare was detected simultaneously both at 
X-ray and TeV energies and the maxima were simultaneous 
within 1 hour, confirming that variations in these two bands are 
closely related (Maraschi et al. 19991. Subsequent more exten- 
sive analyses confirmed these first results also in other sources. 
Note however that the coiTelation seems to be violated in some 
cases, as indicated by the observation of an "orphan" (i.e. not 
accompanied by the corresponding X-ray flare) TeV event in 
the BL Lac IBS 1959+650 (Krawczynski et al. 2004j. 

Although the study of TeV BL Lacs has been the subject 
of great effort in the past (e.g. Bottcher 120041 Georganopoulos 
& Kazanas TUU^ Moderski, Sikora & Blazejowski 200f 
Bottcher & Dermer E002l Tavecchio et al. .2001, Takahashi et 
al.EnOni Kirk, Rieger & Mastichiadis n^^ . a detailed analy- 
sis of the expected correlations of the fluxes in different bands 
is lacking. In this paper we present for the first time a com- 
prehensive study of the time-dependent emission expected in 
the case of the homogeneous SSC model, following the evo- 
lution of electrons and taking into account both radiative and 
adiabatic losses. In particular we focus our attention on the cor- 
relation between the variations observed in the X-rays and in 
the TeV band, which can be quite informative of the processes 
responsible for the observed variations. The present work has 
been stimulated (as described in Sect. 2) by the observation in 
several cases of flares for which a more than linear behavior 
was observed (namely, Fjev °^ -^x-ray' "^^^^ ^ > !)■ first 
sight this does not pose particular problems (an increase of the 
density of the emitting electrons will make the self-Compton 
flux increase more than the synchrotron flux), but we will show 
instead that these more than linear coiTelations are quite dif- 
ficult to reproduce, with standard assumptions, when we con- 
sider that the same type of correlation has been observed not 
only during the rising phases of a flare, but also during the de- 
cay phase (Fossati at al. .2004j . 



In Section 2 we discuss some of the clearest observational 
evidence for the correlation between X-rays and y-rays. In 
Section 3 and 4 we describe the model used to calculate the ex- 
pected correlations in the SSC framework (assuming spherical 
and cylindrical geometries) and we present the results, stress- 
ing the difficulty of obtaining more than linear relations be- 
tween TeV and X-ray variations. In Section 5 we briefly dis- 
cuss the effects that can be induced by considering the finite 
light crossing time in the study of the variabiUty. In Section 6 
we conclude. 

2. Observed correlations 

The correlation between the evolution of the X-ray and the y- 
ray emission of high-energy peaked BL Lac (HBL) objects 
may provide important information on the emission mecha- 
nism of such sources and the physical mechanisms producing 
the variability. However, to investigate such a correlation it is 
necessary to compare two different light curves, obtained in the 
same period of time by two different instruments. The sampling 
rate for both light curves must be similar to provide a reliable 
result. Since the variability timescale in HBL objects is very 
short, the sampling of the light curves should also be dense, 
with several data points per day. Moreover, additional simulta- 
neous spectral observations are required to provide basic infor- 
mation about the evolution of the spectrum. The last condition 
is very important for the modeling of the source activity. For 
example it is important to know if the X-ray synchrotron emis- 
sion (Fs) is generated below or above the vFs{v) peak. It is 
difficult to meet all these conditions. Therefore, there are just 
a few cases where the correlation between the X-ray and the 
y-ray activity can be investigated in detail. In this section we 
present three correlations obtained for two very well known 
HBL objects Mrk 421 and Mrk 501, and we also discuss the 
correlations obtained by other investigators. 

In Figure n we show the activity of Mrk 501 observed in 
1997 April. We show in this figure the light curves obtained 
by the Whipple (E>350 GeV), CGRO-OSSE (50-150 keV), 
and RXTE-ASM (3-20 keV) experiments (Figs [l]-a,b,c re- 
spectively, Catanese et al. I1997I I. We show also the evolution 
of the spectral energy distribution (Fig. QJ-d) observed by the 
BeppoSAX (0. 1-200 keV) and CAT (E>250 GeV) instruments 
during this activity (Plan et al l 19981 Djannati-Atai et al. ll999> . 
For this data we calculate two correlations. The first correlation 
is calculated between the fluxes obtained by the OSSE and the 
Whipple experiments. This correlation gives an almost linear 
result (Fig. ^-e). The precise computation gives in this case 
/^Whipple -^osse"^" ^'''^ ^ reduced;^^^ = 0.67. The second cor- 
relation is calculated between the data from the ASM exper- 
iment and the fluxes obtained by the Whipple telescope. This 
correlation gives a quadratic or even more than quadratic result. 
The precise computation gives Fwhippie ^asm 
ducedx^ = 0.65. Note that to compare two light curves, where 
the observational points are randomly distributed in time, it is 
necessary to select points gathered almost exactly at the same 
time. This may practically mean that we are able to correlate 
only a few points in a light curves where we have even dozens 
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Fig. 1. The activity of Mrk 501 observed in 1997 April. The up- 
per panels show the y-ray and the X-ray light curves obtained 
by the Whipple (a), CGRO-OSSE (b) and RXTE-ASM (c) ex- 
periments (Catanese et al. 1997 ). The shaded areas in these fig- 
ures indicate the observing periods of the BeppoSAX and CAT 
instruments. In panel (d) we show the spectra obtained by the 
above mentioned instruments (Plan et al. 1998, Djannati-Atai 



et al. ll999t . In this panel we show also the energy bands used 
to obtain the light curves. The lower panel (e) shows the two 
correlations made between the X-ray and the TeV gamma-ray 
fluxes. For clarity, we multiply the OSSE data in this panel by 
the factor 0.5. The dashed lines in this panel show a template 
for the linear and the quadratic correlation. 



of points. In the alternate approach we may try to interpolate 
one of the light curves. However, this method requires similar 
sampling of the light curves. Moreover, the sampling should 
be comparable or less than the characteristic variability time 
scale of a source. In the case of the OSSE- Whipple correlation 
we have interpolated the Whipple light curve. The interpola- 
tion method in this particular case has a minimal influence on 
the final result because the OSSE and Whipple observational 
points were gathered almost at the same time. There is only 
one point in the OSSE Hght curve (50549-50550 MJD) which 
has no counterpart in the TeV observations and we have ex- 
cluded this point from our analysis. In the ASM- Whipple cor- 
relation we have interpolated the ASM light curve. The time 
shifts between the observational points in these light curves 
are relatively large in comparison to the previous coiTelation. 
However, each point in the ASM light curve has been calcu- 
lated as an average of many basic measurements made dur- 
ing a period of about one day. Therefore, the interpolation in 
this case also provides meaningful results. Moreover, our re- 
sults are in agreement with the previous correlation made by 
Djannati-Atai et al. il999l l. They obtained a more than lin- 
ear but less than quadratic correlation between the BeppoSAX 
and the CAT telescope fluxes. Note that the energy range of 
the BeppoSAX observations (0. 1-200 keV) covers the spectral 
bands used by the ASM and the OSSE experiments (Fig.[n-d). 
Therefore, the correlation mentioned above gives an average 
of the two correlations presented in this work. The source was 
observed at that time also by other instruments. Krawczynski 
et al. (2000 1 correlated the evolution of the X-ray emission at 
3 keV with the evolution of the y-ray radiation at 2 TeV. The 
correlation they obtained seems to be more than quadratic (see 
Figure 6-a in their paper). This is very similar to the correlation 
presented here between the ASM (2-10 keV) and the Whipple 
observations. The second correlation performed by the team 
mentioned above concerns the evolution of the X-ray emission 
at 25 keV and the evolution of the y-rays at 2 TeV. This corre- 
lation gives an almost quadratic result (see Figure 6-b in their 
article). 

The correlations discussed above show that there was sig- 
nificant difference between the evolution of the synchrotron 
emission before the vF^iv) peak and the evolution around the 
peak. It means that the correlation may depend on the position 
of the spectral bands (below, around or above the vFs/dv) peak) 
used for the calculations. The comparison between our corre- 



lations and the correlation done by Djannati-Atai et al. J1999> 
indicates that the slope of the correlation may also depend on 
the width of the spectral band. 

In Figure |2l we present the activity of Mrk 421 ob- 
served in 2000 February. We show the light curves gathered 
by the HEGRA (E>lTeV), RXTE-PCA (3-20 keV) and the 
RXTE-ASM (2-10 keV) experiments (Figs [T]-a,b,c respec- 
tively, Krawczynski et al. .200 I j . We show also the spectral 
energy distributions observed at that time (Fig. ^-d). We cal- 
culate for this source the correlation between the PCA and 
the HEGRA fluxes. To obtain the coiTelation we have inter- 
polated the HEGRA light curve. However, the interpolation 
as in the case of the OSSE-Whipple coiTelation for Mrk 501, 
has minimal influence on the result because the observations 
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in both light curves were made almost simultaneously. In this 
case we have obtain unclear results; the detailed fitting gives 
■Fhegra -^pcA*" '^ with a poor reduced of 1.7. Note that 
also the correlation between the ASM and the HEGRA data, 
which we do not present here, does not provide an informative 
result. This example shows that the relatively long timescale 
observations with one data point per day may not be sufficient 
to obtain meaningful information on the TeV/X-ray correla- 
tion. 

However, Fossati et al. ( I2004> have recently reported very 
precise observations of this source made in 2001 March. The 
object was observed during a period of one week by the RXTE- 
PCA, Whipple and the HEGRA insti-uments. The rate of the 
measurements was relatively high and several flaring events 
were well observed. The most significant seems to be a flare 
observed on March 19 by the PCA and Whipple instruments. 
The TeV/X-ray correlation for this flare is almost quadratic 
(Fxev « f x^*°-^) for the rising and decaying phase. However, 
the correlation for data obtained by the Whipple and RXTE- 
PCA experiments during this campaign gives Fxev Fi^^'^'^-\ 
This indicates that the longterm correlation may give results 
significantly different from the correlation performed for one 
single flaring event. The reason for the difference may be re- 
lated to the fact that for a long term light curve (e.g. a few days 
or even longer) we may correlate at least a few flaring events. 
If for example the activity is generated in the framework of 
the internal shock mechanism (e.g. Rees 1978 Guetta et al. 
I2004> then the flares may overlap each other in time. If we ob- 
serve a decay phase of a flare and during this decay a new flare 
will start to rise then the observed correlation would be com- 
pletely different even if the base correlation for each flare was 
the same. 

Recently Tanihata et al. (2004 1 published a spectral anal- 
ysis of the observations of Mrk 421 obtained during a 7 day 
campaign in 1998. They found a correlation between the TeV 
flux and the synchrotron peak flux Fjey cc F^-^^^^^. 

In the following we will take particular care in looking for 
robust solutions yielding linear or quadratic correlations be- 
tween X-ray and TeV bands. Although so far there are only 
few cases for which these correlations (especially the quadratic 
one) have been clearly established, our attention is motivated 
by the fact that, as we will discuss below, the observations of 
such a coiTelations could impose strong constraints on the sce- 
nario usually considered for the variability observed in these 
sources. In particular the observation of a quadratic relation 
during the decaying phase is particularly intriguing. In fact, if 
the quadratic increase of the synchrotron and IC emission could 
be easily reproduced, increasing the electron density within the 
source, the decrease is more problematic. Escape of electrons 
from the sources appears unrealistic in the small timescale 
shown by variability. On the other hand radiative cooling 
can affect only the high-energy particles. Indeed, the simple 
quadratic relation predicted for the IC emission is strictly valid 
for IC in the Thomson regime. However, the TeV emission of 
the HBL objects is probably generated in the Klein-Nishina 
regime. In Appendix |A| we estimate the physical parameters 
of a source that could generate TeV emission in the Thomson 
limit. The estimation gives unacceptably large values of the 



Doppler factor (S > 1000, Begelman et al. \T^i . Therefore, 
in all models presented in this work we assume that the TeV 
radiation comes from the IC scattering in the Klein-Nishina 
regime. The decline of the cross-section in the Klein-Nishina 
regime (valid for sources with a plausible value of the phys- 
ical quantities) has the effect of decoupling the (high-energy) 
electrons producing the high-energy photons and those (at low- 
energy) producing the optical-IR synchrotron seed photons 
(e.g. Tavecchio, Maraschi & Ghisellini 1998). In these condi- 
tions, since the radiative cooling only affects the high-energy 
particles, a linear relation between X-rays and TeV should be 
expected, since the seed photons can be considered almost con- 
stant during the decay. A plausible alternative, discussed in de- 
tail below, is to admit adiabatic cooling, affecting all the elec- 
trons, irrespective of their energy. 



3. Spherical homogeneous source 

A homogeneous synchrotron-SSC model is frequently pro- 
posed as a possible explanation for the X-ray and the y-ray 
emission of HBL objects. Such model provide a very good 
and simple explanation for the Spectral Energy Distributions 
(SED) observed in X-rays and the TeV gamma rays (e.g. 
Dermer & SchUckeiser 1993 Bednarek & Protheroe [T997I 
Mastichiadis & Kirk 1 1 9971 Pian et al. II998I Katarzynski et 
al. i2001i) . Therefore, we decided to check if such a model is 
also able to explain quite specific properties of the observed 
variability. In the approach presented in this section we do not 
consider Light Crossing Time Effects (hereafter LCTE) in the 
synchrotron radiation field inside the source nor for the total 
observed emission. This means that the results presented in this 
section are strictly valid only if the physical processes, which 
modify the source brightness, are slower than the source light 
crossing time. This condition may not always be correct in the 
case of blazars. However, it allows us to investigate in detail 
the most important physical processes which may modify the 
emission level of the source. Detailed understanding of these 
processes is very important for any more complex modeling. 

3.1. Description of the model 

We assume for this modeling a spherical homogeneous source, 
which may undergo expansion or compression. As a first step, 
to simplify the model and to provide analytical formulae which 
describe the evolution of such a source, we decided to consider 
only four physical processes which may occur during the evo- 
lution. We calculate increase or decrease of the source volume, 
increase or decay of the magnetic field intensity, variations of 
the particle density and adiabatic heating or cooling of the par- 
ticles. In the model we separate as much as possible the math- 
ematical description for each process. This allows us to inves- 
tigate separately the influence of each process on the source 
emission. In this first approach we do not consider the radia- 
tive cooling of the particles. The impact of this process will be 
discussed in the next part of this work (Subsection l3.7> . 
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The time evolution of the source radius (R) and the mag- 
netic field intensity (B) are assumed to be a power law func- 
tions 



(1) 



where to, Rq and Bq are initial time, magnetic field intensity and 
radius respectively. The indices and m are a free parameters 
in our modeling. 

The initial (f - to) distribution of the electron energy inside 
the source is defined by a broken power law 



Tmin 
„0 



y < r. 



brk 



'''brk ^max 



(2) 



where k2 = ^KTbrk-* "' ' ^ Lorentz factor which is equiv- 
alent to the electron energy, ^^^.j^ describes initial position of 
the break, ni and n2 are spectral indices before and above the 
break respectively. In the model we assume that the dominant 
part of the X-ray and TeV emission is produced by the electrons 
with the Lorentz factors around ybrk- Therefore, the values of 
ymin and y^ax parameters are not important for the results of the 
modeling. For all calculations presented in this section we use 
ymin = 1 andyniax = lOl 

The evolution of the electron energy spectrum is defined by 
a minimum 



N(y, t) = min {A^i (y, f), Nojj, t)} 
of two power law functions 



A^i (y, t) = K, (f)y-'" ; K,{t) ^ ''" ( ^) 



3''<' I In 



Niij, t) = K2{t)y-"^ ■ K2(t) = ^2 ( ^) {'f) 



3'V //„\'o(«2-l) 



(3) 



(4) 



where Ki and K2 describe the evolution of the particle density 
before and after the break iytrk)- The exponent 3rd describes 
the decrease or increase of the particle density. The adiabatic 
heating or cooling of the particles is described by the indices 
'"£i(«i|2 - !)■ Note that if we assume for example adiabatic ex- 
pansion or compression of the source with a constant number of 
the particles, then the parameters r^, rj and should be equal 
(e.g. Kardashev ll962l Longair il992i) . This is in fact the most 
realistic case. However, to investigate separately the influence 
of the above mentioned processes, we decided to produce such 
a detailed parameterization. 

Assuming that the electron spectrum is defined by the min- 
imum of two evolving power law functions we can easily find 
the evolution of the break ybrk(0 = y^^^ (k)lt)'° ■ 

3.1.1. Synchrotron emission 

The synchrotron emission coefficient which describes the radi- 
ation of the electrons is given by the well-known formula oc 
KB^"^^^ (e.g. Rybicki & Lightman [T979t where a = (n - l)/2. 
In our particular case we have to define the evolution of two 
emission coefficients 



j?it) oc Kn2(t)Bit)'-"»'-^'^ <x /:i|2fii|2(y) 



to \3''rf+''<'("i|2-l)+m(Q;i|2 + l) 



where describes emission of the low energy electrons (y < 
ybrk), ]] describes the radiation of the high energy electrons 
(y > ytrk) andBi|2 = Bq'"^'.' 

In our modeling we neglect the electron self-absorption 
process which is important only for the emission at the radio 
frequencies which we do not analyze in this work. Therefore, 
we can approximate the evolution of the intensity of the syn- 
chrotron radiation from the spherical source by 

cc R{t)jf(t) <x Roh\2Bi\2 (y) . (6) 

Finally we can write the evolution of the synchrotron flux mul- 
tiplying the intensity by the source surface 



Flit) (X R^(t)l\{t) oc RlhBii^j] 

ii - 3re - 3rd - rjn^ - 1) - m(ai + 1), 

F^Jt) <x R\t)I^(t) <x Rlk2B2{j) 

S2 - 3re - 3rd - ra{n2 - 1) - m(a2 + 1), 



(7) 
(8) 

(9) 

(10) 



where the indices si and 52 describe the evolution of the syn- 
chrotron emission generated by the low and the high energy 
electrons respectively. Note that there is a difference between 
the evolution of the F] flux and the evolution of the F^ radia- 
tion due to the adiabatic processes (ra(ni|2 - 1)) and the evolu- 
tion of the magnetic field intensity {m{ai\2 + 1)). 

3.1 .2. Self Compton emission 



Tavecchio et al. ( I1998> presented detailed studies of the SSC 
emission generated by the electrons with an energy spectrum 
approximated by a broken power law. The IC spectrum in such 
an approach can be divided into four basic components (F'- ■^). 
However, only two of them are dominant in the IC emission 
generated by HBL objects. The first dominant component (F') 
is generated by the low energy electrons (Ki,Ni) and the first 
part of the synchrotron spectrum (/j,Fj). The second impor- 
tant component (F^) is produced by the high energy electrons 
(K2, N2) and the first part of the synchrotron emission (/] , Fl). 
With these results in mind, we can derive the formulae de- 
scribing the evolution of emission coefficients for the dominant 
components 

jl{t)^Ki{t)ll{t)^Rok]Bi\^jj . (11) 



jl(t) (X K2(t)ll(t) 

I f„ \-f,+6rrf+;v,(n2-I)+ro(ni -l)+m(Q-i +1) 

cx Rok,k2B, (-|^J . (12) 

By analogy to the synchrotron emission we can write the 
evolution of the intensity of the IC emission for both the dis- 

112 112 

cussed cases (f) oc R{t) (f). The evolution of the IC flux is 



Note that in order to reduce the number of equations we use in 



,(5) 



this work the following simplification 



flip means oc ai or 
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given then by 

Flit) <x RHt)llit) cx R^ok^^Bil^j^ (13) 
ci = 4r, - 6rj - 2r«(ni - 1) - m(ai + 1), (14) 

^2(0 (X R\t)f^(t) oc /?4^ifeBi(y) (15) 

C2 = 4rf - 6rrf - rJm - 1) - r„(n2 - 1) - m(ai + 1), (16) 

where ci describes the evolution of the IC radiation in the 
Thomson limit and C2 describes the evolution of the IC emis- 
sion in the Klein-Nishina regime. 

3.1 .3. Components of the synchrotron self Compton 
spectrum: an example 

To show that the approximations used in the derivation of the 
equations describing the SSC emission are correct, in Figure|3l 
we report an example of the evolution of the synchrotron and 
the IC emission generated by the expanding spherical homo- 
geneous source, calculated using our numerical code. The re- 
sults presented in this subsection are correct for a wide range 
of values of the physical parameters usually assumed in order 
to explain the X-ray and TeV emission of the HBL objects. 
However, the results are correct only for case where the domi- 
nant part of the synchrotron and the IC emission (the peaks) is 
generated by electrons with the Lorentz factors close to ybrk 
which practically means n2 > 3. The spectra presented in 
Figure|3lare divided into different components. This test shows 
that indeed only two components are important for the high en- 
ergy IC emission for a large amplitude of the variation. Note 
that the components are calculated numerically, therefore the 
spectra shown are slightly different from the curves obtained 
from the analytical approximations by Tavecchio et al. (1998 1. 
For the numerical calculations we use the SSC emission mech- 
anism described by Katarzynski et al. ( 2001 ). In the test we use 
following parameters 6 = 50, Ro ^ 2 x 10"' [cm], Bo = 0.004 

[G], ^1 = 10^ [cm-3], = 1, ybrk = 1.2 X 10^ r^ax = 10^ 

re = rd - Ta = m - I. Note that in an ideal case where the 
source is perfectly spherical and homogeneous and the mag- 
netic field inside the source is well organized the parameter m 
which describes the evolution of the magnetic field should be 
equal to 2rp (e.g. m = 2 for a linear expansion of the source) in 
order to keep the magnetic flux constant. However, practically 
the occurrence of such ideal sources is very unlikely. Usually 
we would expect a more complex geometry, inhomogeneity 
and a turbulent magnetic field inside the source which may give 
a more complex evolution of the magnetic field. For example 
a slower than expected decrease of the magnetic field during 
the expansion may be caused by the turbulent dynamo effect 
(e.g. Atoyan & Aharonian" 1 999'for discussion of this problem). 
Therefore, we assume the parameter m to be a free parameter 
in our modeling. For the linear expansion we assume that the 
value of this parameter may vary from 1 to 2. The value of this 
parameter has no influence on the results presented in this sub- 
section. However, in other tests presented in the next part of this 
work, we had to use m = 1 in order to explain the quadratic cor- 
relation observed in the decaying phase of the flare. Therefore, 



to be consistent and to provide the possibility of comparison of 
the results, in all numerical calculations presented in this pa- 
per we have selected m - I. Moreover, in this test and in all 
other simulations presented in this work we assume that the to 
parameter is equal to 10 Ro/c. This gives the radius expansion 
velocity equal to 0. Ic for r^, = 1 . To change the source emission 
level at least a few times we simulate the evolution for a period 
of time equal to 3 R[)/c. This set of parameters is quite specific 
because we use a relatively small value of the initial magnetic 
field and a relatively large value of the Doppler factor (see es- 
timations made by Tavecchio et al. 1 19981 or Katarzynski et al. 
12)01 ). We assume a small value of the magnetic field in order 
to neglect radiative cooling, which cannot be described easily 
by the analytical formulae. On the other hand, to reach the ob- 
served level of emission, we have to assume a large value of the 
Doppler factor. The Doppler factor value > 50 has been used 
by some authors (e.g. Krawczynski et al. 120021 Konopelko et. 
al. l2003l l to explain the emission of TeV blazars. We will prove 
in the next part of this work (Subsection 13. 7> that the radia- 
tive cooling is indeed negligible for this set of parameters. In 
this test the radiative cooling is negligible for electrons with 
the energy described by the Lorentz factor around ybik which 
are producing the dominant part of the emission observed as 
the synchrotron and IC peaks. On the other hand the radiative 
cooling is still important for electrons with a Lorentz factor 
around jmax- However, the radiation generated by these par- 
ticles is almost completely negligible due to the steepness of 
their spectrum. 

Moreover, for this test and for other calculations presented 
in this work we use the Hubble constant equal to 65 [km s"^ 
Mpc ' ] and the redshift equal to 0.03. In the above calculations 
and in all other calculations of the IC scattering presented in 
this work, we neglect the possible absorption of the y-ray emis- 
sion due to pair production inside the source. This process can 
be neglected if the Doppler factor of the source is of the order 
of ten or more (e.g. Bednarek & Protheroe ll997> . Moreover, 
to calculate the observed TeV flux, we should also consider 
the absorption of this radiation by interaction with the Infrared 
Intergalactic Background. However, in this particular case, we 
are mainly interested in investigating the evolution of the TeV 
radiation with respect to the evolution of the synchrotron X- 
ray emission. Therefore, absorption of the TeV emission by the 
IR background, which is constant in time, does not modify the 
TeV/X-ray flux correlation. 

3.2. Basic cases 

Correlating the X-ray and TeV emission of the HBL objects we 
can distinguish two general cases. In the first case we compare 
the X-ray emission before the synchrotron peak (Fl) with the 
y-ray radiation above the IC peak [F^, e.g. Mrk 501, Figure^. 
In the second case we correlate the X-ray emission above the 
synchrotron peak (F^) with the y-ray radiation above the IC 
peak (Fp, e.g. Mrk 421, Figure|2j. To simplify the discussion 
presented in this work we hereafter call these correlations C2/S1 
and C2/S2 respectively. We also present in this section the esti- 
mates of other possible correlations (Fl vs Fl and F^ vs Fl). 



Katarzynski et al.: Correlation between the TeV and X-ray emission. 



7 



However, we will not discuss them in detail. The more com- 
plex cases of the TeV/X-ray correlation, where the considered 
bands correspond to the emission around the peaks, discussed 
in subsection l3.5l 

For the tests presented in this section we use the same val- 
ues of the physical parameters as in the calculations presented 
in the previous subsection. To simulate different evolution sce- 
narios we change only four model parameters r^., r^, r^, m set- 
ting or 1 as the value for these parameters. This means that in 
principle we can study only expansion of the source. However, 
as long as we neglect the radiative cooling, the results of such 
tests are opposite to the results which we could obtain from the 
compression of the source. Having four parameters and two 
possible values for these parameters (0 or 1) we have in prin- 
ciple sixteen different combinations (2"*). However, one combi- 
nation is not interesting {re - rd - ra - m - 0) and four other 
combinations are not realistic (all cases where we should cal- 
culate the adiabatic cooling without a decrease of the particle 
density). In Table we show eleven other possible combina- 
tions plus one test where we use m -2. 

We discuss here in detail only a few tests which we consider 
to be the most realistic: 

• In our first case (Table ^-a) we assume the expansion of 
the source volume (r^. = 1) with constant particle density 
and magnetic field intensity (r^ = m - 0). This scenario 
could correspond to a constant injection of particles (e.g. 
from a shock region) which extends the dimension of the 
source instead of increasing the particle density. The syn- 
chrotron radiation in this case is proportional to the source 
volume. Therefore, the synchrotron emission produced by 
the low and the high energy electrons increases with time 
as oc t^'' (ii|2 = 3. Eqs|71 The IC radiation in this 
test is also proportional to the source volume. However, the 
IC emission is also proportional to the intensity of the syn- 
chrotron emission (Eq. |6jl which grows linearly with the 
radius of the source {if oc Rjf oc f'). Therefore, the IC 
radiation increases in time as f^''<-+''' (cip - 4). Finally we 
obtain cip/^ip = 4/3 for all possible correlations. 

• The second case (b) is opposite to the scenario just de- 
scribed. We increase the particle density (rj = 1) keep- 
ing the volume and the magnetic field intensity constant 
{r^ = m - 0). A linear increase of the density provides a 
linear increase of the synchrotron radiation. The IC emis- 
sion is proportional to the density of the particles and to 
the intensity of the synchrotron radiation which is also pro- 
portional to the density. Therefore, the IC emission is pro- 
portional to the square of the particle density. This is well- 
known relation. In this particular case we obtain a quadratic 
correlation in all possible cases (cip/sip - 2/1). 

• In the next case (d) we assume an increase of the vol- 
ume, decrease of the particle density and with adiabatic 
cooUng (rg - rd = Ta = I) with a constant value of 
the magnetic field intensity (m = 0) during the evolution. 
Even if it is rather unlikely that the magnetic field can 
remain constant, we describe this case to provide a sort 
of introduction to our last, more complex test. The syn- 
chrotron emission decreases in this test due to the adiabatic 



Table 1. Estimated correlation between evolution of the syn- 
chrotron emission and the evolution of the IC emission for dif- 
ferent evolutionary scenarios. Columns 2-5 show the parame- 
ters used for the tests. Columns 6-9 show the estimated values 
of the correlation in four possible cases (cip/iip)- 
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losses because the decrease of the density is fully compen- 
sated for the increase in the volume if = r^. This de- 
crease depends on the slope (nip) of the particle spectrum 
(Eq.|4}. The synchrotron emission produced by the low en- 
ergy electrons decreases more slowly (r''"*"'"'', s\ - -1) 
than the emission generated by the high energy electrons 
(^f-r„(n2-\) ^ S2 = -4). The IC emission also decreases in 
this test. In the Thomson limit, the IC flux is proportional 
to the increase in the volume (f^'''"^) and the radius (f''"', 
which comes from the intensity of the synchrotron emis- 
sion l] oc Rj\ oc f') and also to the square of the density 
of the low energy electrons including the adiabatic losses 
{K\ oc ,-6'''-2'-"(«i-i)=-8). This finally gives ci = -4. In the 
Klein-Nishina regime the IC radiation is proportional to the 
evolution of the volume and radius as in the previous case 
(f3r,+r,=4-j fjowcver, this time the scattering is proportional 
to the density of the low energy electrons including the adi- 
abatic losses {K\ oc f-3»',/-''»(ni-i)=-4^ which comes from the 
intensity of synchrotron radiation /]) and is proportional to 
the density of the high energy electrons including the adia- 
batic losses {K2 oc t-^';i-,-,(n2-i)=-7y ^ result, the IC flux 
in the Klein-Nishina regime decreases as (c2 - -7). 
The influence of the adiabatic losses provides in this test 
four different correlations ci/ii = 4, ci/i2 = 1, q/^i - 7, 
C2/S2 - 1-75. 

• In our last test (e) we assume an increase of the source 
radius, the decrease of the particle density with adia- 
batic cooling and also decay of the magnetic field (r,, - 
rd - ra - m - Y). The decrease of the magnetic field 
causes a much faster decrease of the synchrotron emission 
(f-'-„(«ii2-i)-m(tv,|,+i)^ ^ _2_5^ J2 = -7) in comparison to 

the previous case. Moreover, this decrease increases the dif- 
ference between the synchrotron emission produced by low 
and high energy electrons. The IC emission decreases faster 
as well. However, the factor (f-™("i +i)=-i s-j yyjjjcjj marks the 
difference to the evolution discussed in our previous test is 
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the same for the emission in the Thomson and in the Klein- 
Nishina regime. This is a consequence of the fact that the 
y-rays in the Thomson and Klein-Nishina regime are pro- 
duced by the scattering of the synchrotron radiation field 
produced by the low energy electrons (/]). Therefore, the 
decay of the magnetic field does not modify the difference 
between the two scattering regimes. In conclusion, we have 
to add the factor -1.5 (with respect to the previous case) to 
the indices describing the evolution of the IC emission, to 
obtain ci - -5.5 and C2 - -8.5. Also in this case we ob- 
tain four different values for the correlations ci/si - 2.2, 
ci/s2 = 0.786, C2/S1 = 3.4, C2/S2 = 1.214. 

For all estimations presented in this subsection we have as- 
sumed constant values for the parameters «i and «2- However, 
we have checked also that for the most realistic scenarios (e & 
/) and for all combinations of the ni value in the range from 1.5 
to 2.5 and the n2 value in the range from 3 to 8, the correla- 
tion C2/S1 is always more than quadratic while the correlation 
C2/S2 is always less than quadratic. We provide also a simple 
mathematical formulae for these realistic evolutions 

, 1 + 3ni + Ini 1 + 2n\ + no rr \r -v 
C2I s\ - — -^^ '_ — - or 2^ TeV vs soft X-rays 

\+3ni+2n2 I + 2ni +n2 u av 

C2/S2 = — :; ; or TeV vs hard X-rays 



3n2 - 1 



m=l 



2m, 



m=2 



3.3. Problem with "the quadratic decay" 

This concerns the quadratic results of the C2/S2 correlations 
during the decay phase of the flare. Within the scenarios pro- 
posed in the previous subsection there is only one solution 
that may explain exactly the quadratic correlation (see column 
marked C2/i2 in Table QJ. However, this well-known test as- 
sumes only variations of the particle density. Therefore, this 
scenario can be easily used to explain the rising phase of the 
flare, where the activity is generated for instance by the in- 
jection of the particles into a constant volume of the source. 
We cannot easily explain the decrease of the density during the 
decay of the flare. In principle such a decrease could be re- 
lated to the escape of the particles into a region where the mag- 
netic field intensity is significantly less than inside the source. 
However, this process does not guarantee that the efficiency of 
the total IC emission will decrease two times faster than the 
efficiency of the synchrotron radiation. The IC scattering may 
also occur efficiently in the region where the magnetic field is 
significantly less, especially close to the source surface (dis- 
tance < 1.5/?), where the radiation field energy density is on 
average only two times less than on the source surface (Gould 
11979 ^. The decrease of the density could be also related to the 
expansion of the source. However, assuming the expansion we 
have to consider also the influence of some additional physi- 
cal process. We have to consider the adiabatic cooling of the 
particles and probably also the decay of the magnetic field. All 
these processes can destroy the quadratic correlation which is 
related to the variation of the density. Each additional physical 
process which can modify the evolution of the synchrotron and 



the IC emission (Fs\c f*''^) in the same way (f*l'^+-*') can destroy 
this correlation. This is a consequence of the simple fact that if 
cj s - 2 then {c + x)l{s + x) + 2. If the influence of the addi- 
tional process is different for the evolution of the synchrotron 
emission (f'^') and for the evolution of the IC radiation {f^^') 
then we have situations where (c + y)j{s + x) should give 2 if 
we want to keep the quadratic correlation. Therefore, assum- 
ing 5 = -1 and c - -2, which is a simple consequence of the 
decrease of the density, we obtain a simple relation y - 2x 
which should be fulfilled by the additional physical process 
to keep the quadratic correlation. For example the adiabatic 
cooling gives x - -ra(n2 - 1) for the evolution of and 
y - -ra(ni-l)-rQ(n2-l) for the evolution of F^. Therefore, the 
relation y -2x can be fulfilled only if ni = «2 which does not 
provide the correct solution. The efficiency of the synchrotron 
emission depends among other things on the slope of the parti- 
cle spectrum and the intensity of the magnetic field. Therefore, 
the decay of the magnetic field gives x - -m{a2 + 1) for the 
evolution of and y - -m{a\ + \) for the evolution of FI. 
Therefore, the relation can be fulfilled if a\ - 2a2 which also 
does not give the correct solution. Moreover, we have to stress 
that the decay of the magnetic field always results in a much 
faster decay of the second part of the synchrotron radiation than 
the decay of the IC emission in the Klein-Nishina regime and 
the Thomson limit as well. In principle this unwanted effect 
can be compensated by the adiabatic cooling which increases 
the decrease speed of the F'^ flux in comparison to the decrease 
of the F^ flux. However, detailed computations show that the 
adiabatic cooling cannot fuHy compensate the influence of the 
decay of the magnetic field; compare the correlations for the 
cases d (r^ — rd — Ka — l,m — 0), e (Kg — rd — ra — m — I and 
I (fe - rd - ra - l,m - 2) in Tabled 

3.4. Change of the Doppler factor 

It is possible that in TeV sources the Doppler factor 6 changes 
during a single flare. Evidence for this come from the require- 
ment of a large (5-value for the inner parts of the jet, gener- 
ating most of the emission {6 > 10, e.g. Dondi & Ghisellini 
.1995, Tavecchio et al. .19981 Kino et al. l2002l GhiselHni et al. 
I2OO2I Katarzynski et al. l2003l or even 6 ~ 50, Krawczynski 
I2OO2I Konopelko et al. l2003> . while instead at the VLBI scale 
(~sub-pc, including satellite VSOP data) there is no evidence 
of superluminal motion of the jet's components (e.g. Finer et al. 
U999t Edwards & Piner l2002l Finer & Edwards 2004 Giroletti 
et al.,i 2004j . This suggests a strong deceleration of the source. 

To analyze the impact of the change of the Doppler factor 
for the observed TeV/X-ray correlation we assume the follow- 
ing evolution for 6(t) 



(17) 



where d is a free parameter. The Doppler boosting effect may 
change the level of the emission [F(v) oc 6^F'(v')] and the 
observed frequency (v oc dv'). For a power law spectrum 
(F' oc v'"") we then have 



F(v,t)^6(tf^''F'(v,t)oc5ot 



F'(y,t) 



(18) 
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Note that the evolution of the flux at given frequency depends 
not only on the value of the Doppler factor, but also on the spec- 
tral index a. Therefore, the evolution of the synchrotron emis- 

112 

sion can be diff'erent below and above the peak (F^' oc v^^'i-). 
Also the evolution of the IC radiation depends on the spectral 

1 12 — 

index (F^^ oc v^^'ukn^ where q-kn ~ 2Qf2 - ai, see Tavecchio et 
al. [T9^ . 

We can describe the IC and synchrotron flux evolution due 
to the change of the Doppler factor by the y and x coeflicients 
respectively. In this way the coiTelation which includes the 
change of 6 becomes c' /s' = (c + y)/(s + x). The change of 
6 will not modify the original (c/s) correlation if y/x = c/s. 
Values y/x > c/s makes the new correlation c'/s' steeper 
(c' /s' > c/s), and values of y/x < c/s makes the new correla- 
tion c'/s' flatter. This simple estimate shows that the decrease 
of 6 can help to solve the problem of the quadratic decay, but 
only slightly, since the change of the slope of the correlation 
(for actual values of the synchrotron and IC spectral indices) is 
only mild. 

The coefficient y for the evolution of the IC emission in the 
Klein-Nishina regime is given by y2 - -3d - 2da2 + da\. 
The parameter x for the evolution of the synchrotron emis- 
sion below (xi) and above (^2) the peak is defined by xip = 
-2>d - da\\2. Setting a\ - 0.5, 02 - 2, and q-kn = 3.5, as 
in the previous modeling, and assuming d - 1, we obtain 
yi/xi - 1.857 and 3'2/x2 = 1.3. For the most realistic evolution 
scenarios, cases e and I in the Table we obtain a decrease 
of the C2/s\ coiTelation slope ({c2/s\ - 3.4 — > Cj/i'j = 2.5]^ 
and [C2/S1 - 2.5 c'^/s'^ - 2.2]/) and a small increase of the 
C2/S2 correlation slope {[C2/S2 = 1.214 —> c'^/ s'^ - 1.25]^. and 

[C2/i2 = 1^4/4 = 

This simple estimate shows that the impact of the change 
of the Doppler factor is much stronger for the C2/ s\ correlation 
than for the C2/j'2 relation. This process does not modify sig- 
nificantly the slope of the C2/52 correlation, therefore it cannot 
solve the problem of the quadratic decay. 

3.5. Emission around the peaks 

Up to now, we have analyzed the evolution of the synchrotron 
and the IC emission in a case were the observed radiation be- 
longs to the spectrum that can be well approximated by a power 
law function. In this case we can well describe the evolution by 
the simple analytical formulae derived in the previous subsec- 
tions. Now we discuss the case where the emission is observed 
at the peaks of the vFs{v) and vFdv) spectra. In this more com- 
plex case, we use the numerical code to calculate the correla- 
tion. Figurel^shows an example of an evolving SSC spectrum 
during the expansion of homogeneous source. We use for this 
test the same set of physical parameters used for the numer- 
ical modeling presented in subsection 13. II For each emission 
process we selected two spectral bands around the vF{v) peak, 
named A (lO'^-lO"* Hz), B (lO^^-lO''' Hz) and X (lO^^'-lO^^ 
Hz), Y (10^^-10^** Hz). Note that most of the TeV emission is 
usually observed in our X-band. Some observations are still 
possible above the frequency 10^^ Hz (see Figures^d and|2l-d) 
in our Y-band but currently we are not able to obtain observa- 



tions in the whole Y-band. However, the spectral index in the 
Y-band is relatively steep. Therefore, radiation observed in this 
band is dominated by the emission which comes from the low 
frequency part of this band which can be observed. In our test 
we calculate only two coiTelations for above described bands. 

The first correlation (A-Y) is calculated between the evo- 
lution of the synchrotron emission (A) in a transition phase 
(from FI to Fj) and the evolution of the IC radiation (Y) in 
the Klein-Nishina regime (F^). This gives an almost quadratic 
correlation. This particular correlation is a kind of transition be- 
tween the coiTelation C2/S1 (3.4 for this particular test) and the 
correlation C2/S2 (which gives 1.214 for this particular case). 
The second correlation (B-X) is calculated between the evo- 
lution of the second part of the synchrotron radiation (B, F^) 
and the evolution of the IC radiation (X) in a transition phase 
from the emission in the Thomson limit to the radiation in the 
Klein-Nishina regime. The correlation in this particular case is 
almost linear. This is the result of the transition from the corre- 
lation described by ci/s2 (equal to 0.786 in this particular case) 
and the correlation C2/52 (equal 1.214inthis particular case). 

The first test presented above shows that a quadratic corre- 
lation can be obtained in some specific cases when we observe 
the radiation emitted close to the vF{v) peaks. However, we 
stress that this result depends strongly on the spectral bands 
used. This effect is clearly visible in Figure |4] where we show 
the correlation (A'-Y) after a small shift of the A-band (A' 
— > 2 X 10'^-2 X 10'^ Hz). The almost quadratic correlation after 
this shift appears significantly less than quadratic. Moreover, 
to obtain the quadratic correlation we have assumed m = 1. 
This means that the decrease of the magnetic field intensity 
was slower than we could expect from the rule which describes 
conservation of the magnetic flux (m = 2). 

3.6. Evolution of the n2 parameter 

In this subsection we analyse the effects related to the change 
of the index of the high energy electrons spectrum (n2)- The 
main reason for this test is to simulate the activity of Mrk 501 
observed on April 1997. The observations obtained at that time 
by the BeppoSAX satelUte (Fig.[T] Plan et. al. 1998 1 indicate 
that the X-ray emission was almost stable around the energy 
0. 1 keV while the emission around 100 keV was very variable. 

To simulate such evolution we change the value of the n2 
parameter in ten steps, assuming n2 = 2.5 at the beginning 
of the simulation and increasing this value by a factor of 0.25 
after each step. Note that the peak of the synchrotron spectrum, 
at the beginning of the test, is generated by the highest energy 
electrons (ymax = 10^, for this test) because «2 < 3. The other 
model parameters are the same as in the previous modeling, 
except ybik which we changed to 3 x 10^ in order to obtain the 
constant emission around 0. 1 ke V energy. The radiative cooling 
is negligible for the physical parameters used in this simulation. 
In the next subsection we discuss in detail this phenomenon and 
we show that indeed this process is negligible. 

In the test we selected the same spectral bands for the cal- 
culation of TeV/X-ray correlation as in the previous modeling 
and we also discus three different correlations. The first cor- 
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relation (A-Y) is calculated selecting the TeV band decaying 
relatively fast and the part of the X-ray spectrum decreasing 
relatively slowly. At the beginning of the simulation the slope 
of the correlation is almost quadratic and becomes 1 .7 towards 
the end of the simulation. In the second correlation (B-Y) we 
selected the parts of the TeV and X-ray spectra decreasing rel- 
atively fast. At the beginning of the simulation the two fluxes 
vary almost linearly. However, at the end of the simulation the 
correlation slope is close to 0.5. The third correlation (B-X) is 
opposite to the first relation (a fast decay of the X-rays with 
relatively slow decay of the TeV emission) and gives the slope 
close to 0.5 at the beginning of the simulation, changing rapidly 
and reaching the value of 0.13 at the end of the simulation. 

This simulation can reproduce the activity of Mrk 501 ob- 
served in April 1997. The correlation between the Whipple and 
OSSE fluxes (Fig. [0 can be explained as the coiTelation be- 
tween the spectral bands where the TeV and X-ray emission 
evolve relatively fast (e.g. the beginning of our coiTelation B- 
Y). The quadratic correlation between the Whipple and ASM 
fluxes can be explained as the correlation between the part of 
the TeV band which evolves relatively fast and the part of the 
X-ray emission which evolves relatively slowly (e.g. our corre- 
lation A-Y). However, the presented examples depend strongly 
on the position of the peaks and on the selected spectral bands. 

3.7. Impact of the radiative cooling 

In all the scenarios discussed so far we have assumed that the 
radiative cooling of the electrons is negligible. Now we assume 
that radiative cooling is important and we discuss its impact for 
the TeV/X-ray flux correlation. 

To calculate the evolution of the electron energy spectrum 
including the radiative cooling we use the kinetic equation 



dN'iy, t) _ d_ 
dt ~ dy 

where 



{[Ccool(Or^ + Cadm(Or] N\y, t)] , 



(19) 



Ccool(f)= ^— [f/B(0+f/md(f)]: 

3 nieC 



Usit) 



Bjtf 



describes the radiative cooling and C^dm - ^a/t describes the 
adiabatic cooling. For simplicity, we assume that the radiation 
field energy density (f/rad) is equal to the magnetic field en- 
ergy density (Ub)- The initial distribution of the electron energy 
spectrum is given by a continuous broken power law 



No(y,t^ to) = kiy 



1 + 



7 



(20) 



The solution of this kinetic equation is given in Appendix |B] 
However, this solution must be converted to a unit volume 
N{y,t) - N*(y,t),(to/tf'''- to be useful for the calculation of 
the emission coefficient (e.g. Kardashev I 1 962> . Since the evo- 
lution of the electron spectrum is quite complex, we use the 
numerical code to check the influence of the radiative cooling 
on the correlation. 

The radiative cooling may cause a fast decay of the high 
energy electrons. This appears as a cut-off in the high energy 



part of the electron spectrum. In Figure|6lwe show two exam- 
ples for the evolution of the spectrum. For the first assumed 
evolution (Fig. |6l-a) the cooling is almost negligible while in 
the second example (Fig.|6l-b) the cooling significantly modi- 
fies the high energy part of the electron spectrum. In the first 
test we use the same values of the physical parameters which 
we used for numerical calculations presented in the subsections 
I3.1ll3.5l and l3.6l In the second test we use significantly different 
values. We assume a larger initial intensity of the magnetic field 
(Bq = 0.004 —> 0.02 [G]) to increase the radiative cooling rate. 
We assume a smaller value of the Doppler factor (5 = 50 — » 30) 
and smaller initial source volume (Ro^2x 10'^ ^ 1.2 x 10'* 
[cm]) to keep the same level of the emission. Moreover, we also 
modify the initial density of the particles {ki = lO^* — > 3 x 10^ 
[cm-'])- 

Figure |^c shows the evolution of the synchrotron and 
the IC spectrum produced by the electron spectrum given in 
Figure |6l-b. The cut-off' in the electron spectrum produces an 
exponential decay of the synchrotron and the IC emission. 
Figure |6l-c also shows the special bands selected for the cal- 
culation of the correlations. In this test we have introduced one 
more X-ray band (C, 10'^ - 10^" [Hz]) to check the coiTelation 
when using the exponentially decaying part of the synchrotron 
emission. 

The first correlation (A-Y), calculated for this test (Fig. 
|6l-d), is analogous to the A-Y correlation computed neglect- 
ing radiative cooling. This is the correlation between the X-ray 
emission in the transition phase (F' —> F^) and the IC radiation 
in the Klein-Nishina regime. It gives, as in the previous test, an 
almost quadratic result. The second correlation B-X was calcu- 
lated between the second part of the synchrotron emission (F^) 
and the IC radiation in the Klein-Nishina regime (F^, C2I si)- 
In principle this correlation should give a value of the correla- 
tion around 1.2 (see Table [Q. However, the radiative cooling 
significantly modifies the evolution of the X-ray spectrum in 
the B-band. Thus, the X-ray emission decreases faster than in 
the evolution of the power law spectrum discussed in the previ- 
ous subsection (see Figure 0}. Therefore, we obtain an almost 
square root correlation in this case. The last correlation (C-X) 
is the most extreme case where we compare the evolution of 
the X-ray emission in the exponential decay part (C) with the 
evolution of the IC radiation in the Klein-Nishina regime (f ^)- 
The exponential decay causes a very fast decrease of the syn- 
chrotron emission. Since this decrease is more than two times 
faster than the decrease of the IC radiation, the slope of the 
correlation is less than a square root. 

The correlations presented in this subsection depend 
strongly on the position of the spectral bands. Moreover, the ra- 
diative cooling which causes the exponential decay of the spec- 
tra can destroy possible quadratic or linear correlations. This is 
related to the fact that the synchrotron emission responds im- 
mediately to the change of the high energy part of the elec- 
tron spectrum while the relevant IC radiation is generated in 
the Klein-Nishina regime, mainly by the electrons with y less 
than 7bik (see Figure |3l-c-d). Moreover, the radiation field used 
in the scattering is not affected by the radiative cooling. To ob- 
tain almost quadratic result for the A-Y correlation we had to 
assume m—\. 
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4. Cylindrical source - pizza or spaghetti 

So far we investigated the evolution of a spherical source. The 
geometry of the source is not important for the synchrotron 
emission as long as the electron self-absorption is neglected. 
On the other hand the geometry of the source may be very im- 
portant for the IC emission which depends strongly on the radi- 
ation field available for the scattering. If the source is strongly 
extended in one dimension (e.g. x) then only the local radiation 
field, proportional to the other dimensions of the source (y, z), 
is important for the scattering. This effect is important for the 
observed correlation especially if the source expands in only 
one or two dimensions. To investigate this effect we approxi- 
mate the emitting region by a cylindrical geometry, described 
by the radius (R) and the length (L) of the cylinder. We use this 
specific geometry to investigate two opposite cases. In the first 
case we assume L «; 7? a "pizza like" geometry (hereafter we 
call this the pizza case). In the second case we use the opposite 
assumption L » /? a "spaghetti like" geometry (hereafter we 
call this the spaghetti case). 

The main reason for such assumptions is that we are look- 
ing for robust solutions which could explain a more than lin- 
ear (or even quadratic) correlation during the decay of a flare. 
We also investigate when such solutions are independent of 
the specific spectral bands selected for the correlation. Using 
the cylindrical geometry we can simulate the expansion of the 
source in one dimension. This reduces the influence of the adi- 
abatic cooling and the decay of the magnetic field. Moreover, 
with this geometry we can use a lower value of the Doppler 
factor (e.g. 6 ~ 25) and we are still able to avoid the impact of 
the radiative cooling (Bq ~ 0.004 [G]). In addition, for some 
specific viewing angles, we can neglect the light crossing time 
effect (e.g. a pizza-like source viewed face on, or a spaghetti- 
like source viewed at 1 /F, corresponding to a viewing angle of 
90° in the comoving frame). 

We assume that the radius and the length of the cylindrical 
source evolve in time as power law functions 



7?(f) = /?o(y) L(f) = io(y)\ 



(21) 



where Rq and Lq are initial radius and length respectively. 
This parameterization gives a volume which evolves in time 
as V(t) = ttRILq (fo/f)"'^, where E = (2r, + k). 

The initial distribution of the electron energy spectrum is 
assumed to be a broken power law function as in the previ- 
ous model (see Eq.|2ji. The evolution of the spectrum is given 
by a minimum of two power law functions as in the previous 
modeling (see Eq. |3}- However, for this modeUng the defini- 
tion of the power law functions substituted for the minimum 
is different according to the different parameterization of the 
cylindrical geometry. The functions are given by the formulae 



A(«l|2-1) 



A^ii2(r, t) = /:i|2(f)r""^; ^112(0 = Arip {'f)" (^) 



where D - (2rd + Id) describes the evolution of the density and 
A = (2ra + la)/^ describes the adiabatic heating or cooling of 
the particles. 



The evolution of the magnetic field intensity for this model 
is defined in the same way as in the previous modeling (see 
Eq.ID. 

The evolution of the synchrotron emission coefficient is cal- 
culated for the low (7']) and high energy electrons (Jj) sepa- 
rately 

jTit) cc /:i|2(f)B(f)""^^' oc ^ipBip (y ) . (22) 

The evolution of the intensity of the synchrotron emission is 
calculated along the source length 



I'JHt) <^ L(t)fJ\t) <^ Lokip_Bip_{j^ 



(23) 



\ -/,+D+A(ni|2-l)+Mj(Q-ip + l) 
f / 

If we assume that the observer is located outside the source at 
some distance on the symmetry axis of the cylinder then the 
observed flux is given by 

Fl%) cc R\t)lf{t) oc RlLokipBi^ij) (24) 
^i|2 ^ E-D- A(ni|2 - 1) - m(an2 + 1). (25) 

The main difference between the previous modeling and 
the current calculations is in the IC emission. The difference 
comes from the energy density of the synchrotron emission 
iJJrad) which is available for the scattering. 

First we describe the radiation in the pizza case where 
the evolution of the emission coefficient in the Thomson limit 
is defined by 7' (f) oc Ki{t) l\{t). The evolution of the IC 
emission coefficient for the Klein-Nishina regime is given by 
jl it) oc K2{t) ilit). In this particular pizza case we use for the 
calculation of the scattering only local radiation field that is 
proportional to the length of the source (Urad °^ ll <^ y '^)- We 
assume that the contribution of the synchrotron emission from 
the other parts of the source is negligibly small because the ra- 
diation field energy density decreases with the square of the 
distance. According to this assumption the radiation field used 
for the calculation of the IC scattering is isotropic. To obtain 
the evolution of the observed flux we have to multiply the evo- 
lution of the emission coefficient defined above by the source 
volume, which gives 

Flit) ^ RlLlk\Bi(^y\ (26) 
cip = 2re + 21, -2D- 2A(ni - 1) - m(ai + 1), (27) 



for the Thomson limit and 

C2P — 2re + 2le - 2D - A(n\ 
— m{ai + 1), 



1)_A(„2-1) 



(28) 



(29) 



for the Klein-Nishina regime. 

In the spaghetti scenario we use only the local radiation 
field for the calculation of the scattering. However, now the lo- 
cal radiation field is proportional to the source radius (Urad ^ 
j\R). Therefore, the evolution of the emission coefficient in the 
Thomson limit is given by 7' (f) oc Ki(t) R(t) j\it). The evolu- 
tion of this coefficient for the Klein-Nishina regime is defined 
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by jljt) oc K2(t) R(t) jl(t). Finally, the evolution of the ob- 
served fluxes are calculated in the same way as in the pizza 
scenario. For the Thomson limit we have 

Flit) (X RlLok^^Bi , (30) 

CIS ^ 3re + le-2D-2A(ni-l)-m(ai + l), (31) 
and for the Klein-Nishina regime we have 

Flit) cc RlLokik2Bi , (32) 

C2S = 3re + le-2D- A(ni - 1) - A(n2 - 1) 

- miai + 1). (33) 

4.1. Specific solutions 

With respect to the spherical geometry, the cylinder is 
parametrized by two more free parameters (Lq, le) which give 
more possible scenarios for the time evolution of the flux. In 
this subsection we discuss only a few possible cases that give 
directly a quadratic correlation between the second part of the 
synchrotron emission (F^, ^2) and the IC radiation in the Klein- 
Nishina regime (F^, C2) during the decay of a flare. 

• Consider a pizza-like geometry and only linear increase 
of the source length He - \, r^, - rd - - Id - h - 
m - 0). The increase causes a linear growth of the syn- 
chrotron emission (F''^ oc t^',s\\2 - 1) and a quadratic 
(c(i|2)P = 2) increase of IC radiation which is proportional 
to the increase of the volume and the radiation field as well 
iUrad jsL oc t''). Thus in all possible cases we can ob- 
tain the quadratic correlation (c(i|2)/>/ii|2 = 2). The same 
test for the spaghetti case gives C(ij2)p/ii|2 = 1 because the 
local radiation field is proportional to the radius and is con- 
stant. Note that the similar test for the spherical geometry 
(r^ = l,rrf = = m = 0) gives cip/^ip = 1.333. 

• In the next test we assume a pizza-like geometry with 
expansion of the radius, decrease of the density and the 
adiabatic losses that correspond to the change of the ra- 
dius (r^ - fd - ra - 1). The other parameters are as- 
sumed to be constant (Z^ - Id - h - m - 0). For 
re - rd the expansion of the volume is compensated en- 
tirely by the decrease of the density. Therefore, the syn- 
chrotron emission decreases only due to the adiabatic cool- 
ing (F? oc 52 = 2.67). The IC emission in the 
Klein-Nishina regime depends on the expanding volume 
it^'''"^), the square of the decreasing density (r'*'^''""'*) and 
the adiabatic losses of the low and high energy electrons 
(f-2r„(«ii2-i)/3=-(0.67|2.67)-) combination of the above pro- 
cesses gives a decrease F^ oc r^^^, which finally provides 
a quadratic correlation ic2p/s2 - 2). However, the corre- 
lation depends on the value of the ni and n2 parameter 
(ni|2 - 2|5, in the above test). We have checked that for 
a given value of the parameter ni in the range from 1.5 to 
2.5 there is always only one value of the index «2 which 
provides the quadratic result (e.g. ni - 2.5 and n2 - 5.5). 
Other correlations for this case, which we do not discuss 
in detail, are described by cipjsx - 5, cip/s2 - 1.25, 
C2p/Si = 8. 



• The quadratic correlation can be obtained also in a three 
dimensional expansion of a pizza-like source (r^ - rd - 
r„ - le - Id - la - However, to get the quadratic 
result it is necessary to assume a constant magnetic field 
(ot = 0) during the expansion and specific values of the 
parameters rti and n2 (e.g. n\ - 2 and n2 - 4). Note that 
similar three dimensional expansion of the spherical source 
(r^. = 1 = = = l,m = 0) gives C2/.S2 = 1-75. 

• The spaghetti scenario also provides a quadratic correlation 
if we assume the expansion of the radius and the length of 
the source (r^ - rd - rg - U - Id - h - However, also 
in this case we have to assume a constant magnetic field 
im - 0) and specific slopes of the electron spectrum (e.g. 
Hi -2 and «2 = 4) to get the right correlation. 

We conclude that the cylindrical geometry can provide 
some quadratic solutions if we can assume a constant or al- 
most constant value of the magnetic field during the evolution 
of the source. Moreover, we have also to assume that the im- 
pact of the radiative cooling is negligible. Otherwise with this 
particular geometry the explanation of the quadratic result of 
the correlation C2/S2 is rather problematic. 

5. Light crossing time effects (LCTE) 

We distinguish two different types of LCTE. The first effect, 
which we call the internal LCTE, is important for the calcula- 
tion of the evolution of the radiation field inside the source. 
This effect is especially important for precise calculation of 
the IC emission. The second eff'ect, which we call the exter- 
nal LCTE, is important for the observer. We investigate in de- 
tail in this subsection this second effect, which seems to have a 
stronger impact for the correlation between the TeV and X-ray 
emission. Note that recently Sokolov, Marscher & McHardy 
(2004 1 proposed a new model to explain the rapid multifre- 
quency variability of blazars. They investigated in detail the 
influence of LCTE for the observed emission. This model pro- 
vide especially precise description of the IC scattering in a case 
where the internal LCTE is important. 

To check the influence of the external LCTE for the cor- 
relation we have used the model proposed by Chiaberge & 
Ghisellini ( 1999 1. They assumed that the source is created by a 
shock wave which accelerates the electrons. The electrons are 
injected from the shock region into a volume of the jet where 
they generate most of the synchrotron and the IC radiation. The 
geometry of the emitting region is approximated to be cubic 
and the region has been divided into small homogeneous cells. 
In figure we sketch the evolution of such a source. On the 
right we show the intrinsic evolution. On the left we show the 
evolution of the part of the source that can be observed in the 
comoving frame at 90° with respect to the shock velocity. The 
number in each cell indicates the age of the particles inside the 
cell. 

If the external LCTE is not important then the observer can 
see immediately the changes of the total source structure. This 
condition is true for this model if the shock velocity (ySj) is 
significantly less than the speed of fight. Otherwise only part 
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of the source volume can be observed by the observer at given 
time. 

We examine two theoretical scenarios of the source evolu- 
tion. In the first we assume that the emission of a single cell is 
constant in time. For the second scenario, opposite to the pre- 
vious one, we assume a very fast flux decay within a single 
ceU. 

5.1. Constant emission of a single cell 

The assumption of constant emission requires the same age of 
the particles in all cells. This means that in cell presented in 
Figure we should have the same number-1. To check the 
influence of the process for the evolution of the emission we 
compare the intrinsic evolution of the source with the observed 
evolution. 

For the intrinsic evolution of the source the constant injec- 
tion of the cells gives a linear increase of the synchrotron emis- 
sion, which is proportional to the increase of the source vol- 
ume. If we neglect the internal LCTE, the IC radiation also will 
increase proportionally to the increase of the volume. However, 
this assumption is not precise and we will discuss the problem 
of the internal LCTE at the end of this section. Immediately 
after the the injection stops (for example due to destruction of 
the shock wave) we should observe a constant emission of the 
source where the absolute level of the emission is proportional 
to the created volume. 

If we consider the external LCTE then observer at the be- 
ginning of the evolution (f - At) will receive only the emission 
from the cell which is nearest to him (bottom row of the cells). 
At a later time the observer will start to receive successively 
emission produced by the cells located at larger distances. If 
the emission of a single cell is constant in time (number 1 in all 
observed cells) then the observed emission increases as a sum 

t^-'^' (34) 

k=i ^ 

where n is the number of cells in the bottom row. Note that this 
formula is correct only for the source which at the end of the 
injection is a square (nx« cells) and is correct only for the injec- 
tion phase. This simple estimation indicates that the observed 
flux will increase in time as if « » 10. This result is signifi- 
cantly diff'erent from the intrinsic evolution where the increase 
of the emission in the injection phase was linear. Therefore, 
this simple test shows that the external LCTE may significantly 
change the evolution of the observed emission. Note that in this 
test the level of the emission increases also after the end of the 
injection. However, in this phase the increase is slower than 
and after some time (2 x Z/c where / is the dimension of the 
square source) the flux level becomes constant, and equal to 
the level of the intrinsic emission after the end of the injection. 

5.2. Very tast decay of a single cell 

We now assume that the flux, within a single cell, decays very 
fast. The "life time" of a single cell is At. 



If the external LCTE is not important we observe constant 
emission of the source as long as the shock is injecting parti- 
cles. The total volume of the source during the whole injection 
process is represented only by the cells next to the shock front 
(e.g. Figure0-d). The source disappears abruptly after the end 
of the injection. 

When the external LCTE is important, at the beginning of 
the injection the observer will see only the radiation produced 
by the bottom row, reduced in this particular case to a single 
cell. At a later time the radiation produced by the cells at larger 
distances will arrive at the observer. This means that the ob- 
server will see a linear increase of the emission. The duration 
of this increase depends on the number of cells {k) in the col- 
umn next to the shock front. After the time kx At the observed 
emission should be constant. 

After the end of the injection we should observe a linear de- 
crease of the emission with a duration equal to the light cross- 
ing time from the back to the near part of the source. This is the 
opposite behavior to the linear increase observed at the begin- 
ning of the simulation. 

There are two important conclusions from this test. The first 
is that even if the intrinsic emission of the source is constant 
for some time at the beginning of the evolution, the observer 
will see a linear increase of the radiation. The second conclu- 
sion concerns the decay of the source. Even if the source "dies" 
abruptly in a very short time, the observer will see a linear de- 
cay of the emission for a time corresponding to the dimension 
of the source. The conclusions obtained in this and in the pre- 
vious subsection are strictly valid for the specific adopted ge- 
ometry. In the next subsection we discuss more possible impli- 
cations of LCTE. 

5.3. LCTE - general conclusions 

The model used to explore the influence of LCTE for the evo- 
lution of the source emission gives the possibility to investi- 
gate only the external LCTE. However, the conclusion obtained 
from the simple test which we have performed with this model 
seems to be quite strong. If the emission produced by the total 
volume of the source evolves in time as t'' then the observed 
emission evolves as f ''^'. Besides being confirmed by our sim- 
ple analytical tests, this very simple rule is also confirmed by 
numerical calculations assuming several different decay rates 
of a single cell emission. However, this result is strictly valid 
only for the proposed geometry, source evolution and specific 
position of the observer (at 90°). The model used cannot pre- 
cisely describe the evolution of the radiation field inside the 
source (internal LCTE) and cannot describe the emission of 
the expanding source. This would require a more sophisticated 
numerical modeling. 

If we assume that the external LCTE indeed increases by 
one the index which describes the evolution of the emission, 
then we can analyze the impact of this process for the ob- 
served congelation. If for example the correlation between the 
synchrotron and the IC emission is intrinsically quadratic then 
the observed correlation will not be quadratic. This is a sim- 
ple consequence of the rule already discussed in subsection 
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lO (if cIs = 2 then (c + x)l{s + x) 2). On the other 
hand if we observe a quadratic correlation it means that the 
intrinsic change of the IC emission must be much larger that 
the change of the synchrotron radiation. For example we need 
Fs f and f '"^ intrinsically to obtain a quadratic corre- 
lation after the basic transformation due to the external LCTE 
((c = 3 + l)/(s = 1 + 1) = 2). The transformation indicates 
that only some specific intrinsic evolutions of the synchrotron 
and the IC emission (e.g. c - l&s-O, c-l>&.s — 1, 
c = 5 & s = 2) may provide a quadratic correlation. Similar 
arguments apply to any observed coiTelation slope. 

The significantly faster evolution of the intrinsic IC flux 
with respect to the evolution of the intrinsic synchrotron flux 
seems difficult to explain. If we neglect the internal LCTE, 
then for the calculation of the IC scattering inside a single cell 
we use only the synchrotron radiation field produced within 
this cell. With this assumption the evolution of the IC and syn- 
chrotron radiation will be exactly the same (c - s). However 
this approach is of course not always exact, since electrons in- 
side a given cell may scatter the seed photons produced also by 
the surrounding cells. The cell at the center of the source is the 
one most affected by this effect. For such a cell we could have a 
linear increase of the radiation field during the injection phase. 
The cells affected the least will be the ones at the corners of the 
source (for them the effect should be four times less than for 
the cell in the center). Therefore, the evolution of the intrinsic 
IC emission of the whole source should be more than linear but 
less than quadratic (Fc oc f' ^ ^) which after the transforma- 
tion gives 2 < c < 3. This estimation requires constant intrinsic 
synchrotron emission to explain a quadratic correlation in the 
observer frame. 

To conclude: i) if the intrinsic evolution of the source pro- 
duces a given correlation between the synchrotron and IC emis- 
sion, then the external LCTE yields a correlation with a differ- 
ent slope; ii) if we observe a given coiTelation and the external 
LCTE is important, then the intrinsic relative IC/synchrotron 
change must be much stronger than what is observed. 

6. Summary and conclusions 

We have presented a detailed study of the expected correlation 
between the variations observed in the X-ray and TeV bands 
in HBL, in the context of the widely-used homogeneous SSC 
model. This work has been stimulated by the observations of 
both linear and quadratic correlations in the few cases for which 
the available data are suitable for a detailed study of the corre- 
lation in single flaring events (Section|2j. 

First we addressed the problem in the context of the widely 
used spherical geometry (Section 0, presenting analytical re- 
lations valid when the radiative cooling can be neglected and 
numerical results valid in general. While a quadratic correla- 
tion during the increasing phase of a flare is easy to reproduce 
through an increase in the density (due, for instance, to a con- 
tinuous injection of new particles in the emitting volume), we 
found that the same correlation observed during the decreas- 
ing phase poses a difficult problem. A case close to a quadratic 
relation can be obtained only in the rather physically implausi- 
ble case of adiabatic expansion and a constant magnetic field. 



When the observational bands comprise the peaks of the SEDs 
the situation is more complex even if a close-to-quadratic rela- 
tion can be obtained, the solution strongly depends on the exact 
position of the bands and a small change in the observational 
limits inevitably would change the coiTelation. 

The main conclusion that we can derive from the first part 
of our study is that in order to to get a quadratic (or, even more 
general, more-than-linear) correlation between X-rays and TeV 
we need rather special conditions and/or fine-tuning in the tem- 
poral evolution of the physical parameters. For special choices 
of the spectral bands under study and of the parameters de- 
scribing the evolution we can get even more-than-quadratic re- 
lations. In all the cases these solutions are quite "delicate" a 
small change in the observational band and/or in the parame- 
ters inevitably changes the correlation. 

Next (Section |4} we investigated the changes related to a 
source geometry, assuming that the emitting region is a cylin- 
der. Synchrotron emission is not affected by the actual shape 
of the emission region (as long as absorption effects within the 
source are negligible). On the other hand the value of the ra- 
diation energy density strongly depends on the geometry, af- 
fecting the SSC emission. We distinguish between two cases 
(called pizza or spaghetti) depending on the ratio between the 
length and radius of the cylinder {R/L>loiR/L<l respec- 
tively). We found that, similar to the spherical case, quadratic 
solution can be found only assuming that the magnetic field is 
not affected by the expansion. Moreover, radiative cooling can 
destroy the correlation. 

Although light crossing time effects have not been directly 
taken into account in the calculations, we have briefly dis- 
cussed possible effects of this phenomena (Section|5|i. We show 
that the external LCTE can "weaken" the intrinsic (within the 
source) correlation. On the other hand if we observe for exam- 
ple a quadratic correlation and the external LCTE is important 
then the intrinsic IC/synchrotron change must be stronger than 
quadratic (e.g. c/s = 3/1). 

As a last step we discussed the possibility that the SSC 
occurs in the Thomson regime. This would easily produce a 
quadratic correlation between X-rays and TeV, since the seed 
photons for the IC scattering are produced by the same elec- 
trons. However this would require rather implausible condi- 
tions for the emission region (extreme Doppler factors, 6 ~ 
1000, and an extremely compact source, see AppendixlAl. 

The general conclusion is that if the quadratic correlation 
during single flares will be found to be common, the sim- 
ple homogeneous SSC model will face a severe problem. As 
stressed many times, the main difficulty is to explain the de- 
caying phase, therefore acceleration processes and/or particular 
injection mechanisms cannot contribute to solving the problem. 
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Appendix A: Large Doppler factor - the Thomson 
limit 

The main difficulty in explaining a more than linear correla- 
tion between the TeV and the X-ray fluxes is associated with 
the fact that in TeV sources the seed photons to be scattered 
at high energies are those below the peak of the synchrotron 
emission, because the Klein-Nishina decline of the scattering 
cross section with energy implies a small efliciency for the scat- 
tering between photons at the synchrotron peak and electrons 
with random Lorentz factor close to ybrk- This "decoupling" 
between scattering electrons and scattered seed photons intro- 
duces the problems we have discussed. 
Indeed, Tavecchio et al. 



where 



I1998t have demonstrated that in 
the single zone, homogeneous SSC model, the knowledge of 
the peak frequencies Vj and v^, the flux at these frequencies, 
and the variability timescale is sufficient to determine the phys- 
ical parameters of the source. In this case one associates the 
variability timescales fvar with the size R of the source, through 
R = cfvar^- Application to Mrk 421 and Mrk 501 using this pre- 
scription results in the finding that both sources are suffering 
from Klein-Nishina effects. 

We now relax the assumption R = c?var<5, and require in- 
stead that the scattering of X-ray photons with the electrons 
producing these same photons (by the synchrotron process) oc- 
curs in the Thomson limit. This will have the advantage that 
any variations of the number of electrons emitting in the X-ray 
band will automatically produce a quadratic variation in the 
TeV band. 

To this aim we define as the random Lorentz factor of 
those electrons producing synchrotron photons of frequency Vx, 
and require (Thomson limit) 



hv 



■_y± , 

2 5 



1 



Since Vj ^ 3.6 x 10*Byj-5, the above condition yields 

B6^ > 



(A.l) 



(A.2) 



3.6 X 10^ 

If the scattering occurs mainly in the Thomson limit, then the 
ratio between the self Compton and the synchrotron peak fre- 
quencies is Vc/vs = (4/3)7^^1^. Since the synchrotron peak fre- 
quency is given by v., ^ 3.6 x IQ^ By^^^d, we have 

„2 



B6 = 



(A.3) 



2.8 X lO^Vc 

The third condition concerns the synchrotron and self Compton 
peak fluxes. The intrinsic (comoving) synchrotron radiation en- 
ergy density is given by = L^yn/A-nR^cd'^. Since, in the 
Thomson regime, we have L^/Ls = U'^y^-^/Us, we have the third 
requirement 

^l/2 



Bd'^R = 



cLc 



(A.4) 



Eauations lA.2llA3l and lA.4l form a closed system of equations 
in the unknowns B, 6 and R. The solutions are 



B 



R 



CLr 



1/2 



1 



(A.5) 



h 



m.c^ 3.6x10^' 



b = 



1/2 



2.8 X iCv, 



(A.6) 



For illustration, consider the case Vj = 1 keV, Vc - 0.5 TeV, 
Vx - 10 keV, and equal synchrotron and self Compton lumi- 
nosities, Lsy„ - Lc - 10"^^ erg s"'. Then we obtain B = 0.14 G, 
6 - 1200 and R - l.2x 10'^ cm. With these parameters, the 
cooling frequency after one hour (observed time) is of the or- 
der of 5 X lO"" Hz. The minimum variability timescale, fixed by 
the light crossing time, is fyai- ~ R/c6 - 0.03 s. As we can see, 
the required beaming factor is rather extreme for this reason we 
consider this solution very unlikely. 

Appendix B: Solution of the kinetic equation 

Analytical solution of the kinetic equation used in subsection 
I3.7l is given by following formula 



N:(y,t) = ki{Si(y,t)/yrS2(y,t) 

Si(y,t) = [/i(f)-r/2(f)], 
SiiyJ) = exp 



1 + 



1 



(B.i) 



r {rCadia(f')/2(f) - rCadia(f')/2(f') 
Cadia(f')/l(0-2r/l(f')Ccool(f')] 

^h(t')-yh(i) + h{i)\At' , 



l\{x) - exp 



Cadia(3')d3' 



h{x) 



Ccooi(3')/i(3')d3', 



This solution describes evolution of the electron energy spec- 
trum only for the case where the initial spectrum is approxi- 
mated by the continuous broken power law function (Eq.l20>. 
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Fig. 2. The activity of Mrk 421 observed on 2000 February. 
The upper panels show the y-ray and the X-ray light 
curves obtained by the HEGRA (a), RXTE-PCA (b) and 
RXTE-ASM (c) experiments (Ki-awczynski et al. IMJU . 
The shaded areas in these panels indicate the observ- 
ing periods for the spectra obtained by the RXTE and 
HEGRA instruments presented in panel (d). The lower 
panel (e) shows the correlation between the X-ray and 
the TeV y-ray fluxes. The dashed lines in this panel show 
a template for the linear and the quadratic correlation. 
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Fig. 4. The upper panel shows the evolution of the SSC emis- 
sion of an expanding spherical homogeneous source. The areas 
indicated by the capital letters show the spectral bands selected 
for the calculation of the correlation. The lower panel shows 
the two main correlations (bold solid lines with arrows which 
indicate the direction of the evolution) calculated for the evolu- 
tions presented in the upper panel. The bold dashed line in the 
lower panel (A'-Y) shows the correlation obtained after a small 
shift of the A-band. Thin dashed lines show a template for the 
correlations. 
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Fig. 3. The evolution of the synchrotron and the IC emission of 
an expanding spherical homogeneous source. Panel (a) shows 
the evolution of the electron energy spectrum. The spectrum 
gives the synchrotron radiation shown in panel (b). The evolu- 
tion of the IC spectrum is shown in panels (c) and (d). To sim- 
plify the figures we show only the initial and the final spectra. 
Note that the initial spectra in all figures are shown by bold- 
solid lines while the spectra generated at the end of the simu- 
lation are shown by bold-dashed lines. The basic components 
of the emission are shown by various thin lines. A detailed de- 
scription of the components and the parameters used for the 
test are given in subsection l3.ll 
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Fig. 5. The upper panel shows the evolution of the SSC emis- 
sion of the source where only the slope «2 of the high energy 
part of the electron spectrum was modified during the simula- 
tion. To calculate the correlations presented in the lower panel 
we selected the same spectral bands as in the previous mod- 
eling. Thin lines in the lower panel show a template for the 
correlations. 



cn 
o 



E 

o 

cn 

0) 



CT> 
O 




20 22 24 
109,0(1^) [Hz] 



o 
cn 

<D 



cn 
o 




-13 



-15 -14 -13 -12 -11 
logio(^s'(f)) [erg 



-10 



Fig. 6. The impact of the radiative cooling for the coiTelation. 
Panel (a) shows the evolution of the high energy part of the 
electron energy spectrum for the physical parameters which we 
have selected for the modeling in the subsection l3.ll Panel (b) 
shows the evolution of the electrons spectrum for the parame- 
ters used in the subsection l3.7l The evolution of the synchrotron 
and the IC emission is shown in panel (c). We also show the 
spectral bands selected for the calculation of the correlations. 
Panel (d) shows three different correlations calculated for the 
tested evolution. 
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Fig. 7. The evolution of the source according to the model pro- 
posed by Chiaberge & Ghisellini (1999). The right column 
shows evolution of the total volume of the source. The left col- 
umn shows the cells that contribute to the source emission ob- 
served at given time by the observer located at 90° (bottom of 
the figure) with respect to the shock velocity (J3s) in the co- 
moving frame. The number in each cell indicates the age of 
the particles. Note that at time 4Af the injection process was 
stopped. 



